Instability statistics and mixing rates 
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We claim that looking at probability distributions of finite time largest Lyapunov exponents, and 
more precisely studying their large deviation properties, yields an extremely powerful technique to 
get quantitative estimates of polynomial decay rates of time correlations and Poincare recurrences 
in the -quite delicate- case of dynamical systems with weak chaotic properties. 
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I. INTRODUCTION 



A general problem of dynamical systems theory con- 
cerns both the evaluation of time or space averages, and 
an understanding of how finite order estimates of such 
averages converge to their asymptotic limit. The second 
feature is tightly connected to a quantitative estimation 
of mixing rates, that is how time correlations asymptoti- 
cally decay. Though in extremely simplified models (sim- 
ple Markov chains, lattice models with a finite transition 
matrix) the connection between convergence of finite or- 
der estimates and mixing properties may be established 
quite easily, the issue in more general contexts is much 
more complicated, and many facets of the problem still 
remain as open problems. In particular when weakly 
chaotic systems are considered it is expected that stick- 
ing to regular structures in the phase space severely de- 
grades mixing properties, and critical, polynomial decay 
of temporal correlations can be observed. We remark 
that mixed systems are thought to be generic and 
that slow polynomial decay may influence deeply deter- 
ministic transport properties as the Kubo formula 
suggests. We here propose to look at large deviation 
properties of finite time estimates as an efficient tool to 
get quantitative information about mixing properties of 
the system. We emphasize that the idea of a relationship 
between the tails of finite time distributions and mem- 
ory effects is not new: for instance inspection of the tails 
of finite time Lyapunov exponents was used to get in- 
formations about qualitative changes in the dynamics of 
coupled maps in Q, the detection of small regular is- 
lands in 0, 01 , or motion of a particle in a random time- 
dependent potential [1] (see also for an example in 
hyperbolic dynamics). A rigorous analysis, inspiringthe 
present work, was presented in [1,11: ori ginally [8| a 
class of onc-dimensional maps / was taken into account 
(the theoretical aiialysis was extended to higher dimen- 
sional systems in [lO[ : if we denote by A the Lyapunov 
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exponent, and by P,i(A„) the distribution of finite time 
estimates 
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then, by fixing a threshold A such that < A < A we 
may estimate the fraction of initial conditions yielding 
an estimate below the threshold 



dXn Pn (An). 



(2) 



The quantity A^^(n) asymptotically vanishes if the sys- 
tem is ergodic (as we always suppose): for weakly chaotic 
system it may however decay polynomially 
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in such a case a bound is proven [8| for correlation decay 

I 



C{n) < 



(4) 



where, as usual, 
Cin) 
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where and tjj are taken in a suitable class of smooth 
observables and fi is the invariant measure of the sys- 
tem. Actually in Q A4j^{n) is defined in a slightly differ- 
ent way (a full equivalence is established only via further 
assumptions) : our definition accomplishes a twofold pur- 
pose: on the one side it allows a comparison with large 
deviation estimates in Q , and on the other side it yields 
an easily computable quantity. 

This actually leads to the main point of our paper: we 
argue that scrutinizing the way in which Mj^{n) decays 
provides an extremely efficient way of studying quantita- 
tively the decay of correlations, a key issue in the anal- 
ysis of fully and weakly chaotic dynamical systems. As 
a matter of fact direct quantitative estimates of mixing 
rates are known to be quite delicate numerically |12| , and 
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many efforts have been devoted to devise how to tackle 
the problem by alternative approaches, like return time 
statistics [H, [13, d , which can be rigorously shown 
to yield the correct answer in Id intermittency \W^ . 

The plan of the paper is as follows: firstly we suggest, 
in view of recent results 0, that the estimate (|4|) is not 
optimal, and the decay properties of A4^{n) and C{n) 
arc the same; then we check such a conjecture for a class 
of one-dimensional intermittent maps, where polynomial 
decay rates of correlations are known exactly. We then 
consider an intermittent area preserving map, which is a 
prototype example of intermittency in higher dimensions. 
We finally reinvestigate, by this technique, the problem 
of generic correlation decay in maps with a mixed phase 
space: in all cases the results are coherent with exact 
results, and corroborate proposed conjectures, in a clean, 
controlled way. 



II. LARGE DEVIATIONS AND CORRELATION 
FUNCTION 



The way in which we defined Mj^{n) ^ is in the form 
of a large deviation result (T^]: a remarkable connec- 
tion between large deviations properties and correlation 
decay was proven in (see also [3): if we consider 
a system for which correlation of smooth functions de- 
cay poly nomially (the smoothness is essential see 
also [i^) C{n) ^ then the following results hold, for 
(Cesaro) finite time averages of an observable (f>: 



fl X 



-15]0(/W(x-)) 



fc=0 
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(5) 



/i being again the invariant measure, while 4> denotes the 
phase average of the observable. In one dimension, by 
taking (/)(a;) = In \ f'{x)\ this suggests that the bound ^ 
could be sharpened, in such a way that correlation decay 
and large deviations are characterized by the same power- 
law exponent: on the other side in dimensions higher 
than one finite estimates of Lyapunov exponents (com- 
puted for instance by selecting the largest eigenvalue of 
the product of jacobians along a trajectory) lack composi- 
tion structure of Cesaro sums. Our claim is that, despite 
these provisos, the asymptotic decay of M.y^[n) quantita- 
tively reproduces exactly the decay of correlations, and 
inspection of such a quantity provides an excellent tool 
to investigate the weak-chaos regime. 

Before applying such a method to a number of dy- 
namical systems, we point out that the way in which we 
consider the exponent ^ is as the lowest possible value 
in a suitable class of (zero mean) smooth functions: it 
is well known that specific choices of observables may 
yield faster decay (see [20| for a rigorous example in one- 
dimensional intermittent dynamics). 



A. One-dimensional system: the Pikovisky map 

Our first benchmark tool is represented by Pikovsky 
map [2ll | Tz^ which is implicitly defined by 
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Yz 



'l + T,{x)Y 



Q <x < l/(2z), 



T,{x) + —[l~T,{x)Y, l/(2z)<x<l; 



(6) 



while for negative values of x, the map is defined as 
Tz{~x) ~ —Tz{x) [2l|. Such a map has remarkable fea- 
tures, as in X = ±1 two marginal fixed points of the 
Pomeau-Manneville type are present (and z is the corre- 
sponding intermittency parameter), while for any value 
of z the invariant measure is Lebesgue (2]| . as it can be 
seen by direct inspection of the Perron-Frobenius oper- 
ator. For this map the correlations decay polynomially, 
with an exponent that depends upon the intermittency 
parameter z 



C{n) 



,-1/(^-1) 



(7) 



This decay is very well reproduced by numerical simula- 
tions for M^{n), as we see from fig. ([Ij. We remark that 
the numerical simulations for the Pikovsky map were per- 
formed by using 10^ initial conditions uniformly chosen 
in the phase space [—1, 1]. 
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FIG. 1: Numerical data (symbols), and regression fits (lines) 
of power-law decay of M^{n) for Pikovsky map for three dif- 
ferent values of the intermittency exponent. The fit for the 
case z = 1.5 starts at n = 3000. Theoretical values for ^, from 
(III), are ^ = 2, ,t = 1 and ^ = 1/3. 

When we want to check how well correlation decay is 
reproduced by large deviations, expressed by the decay of 
Ai)^{n), we have to be careful to take the large deviation 
parameter A sufficiently far from the centre of the distri- 
bution: even if the general expectation is the asymptotic 
decay law will not depend upon the cutoff, finite time es- 
timates require that A4 picks up only contributions in the 
tail of the distribution. As a matter of fact any choice of 
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the cutoff A will determine a time scale such that only 
after such a transient the asymptotic decay is reproduced 
correctly. This is well illustrated in fig. ^ where were 
chosen three values of A and we obtained the asymptotic 
decay for all cases, but with different time scale r^. 



asymptotic phase average, from a physical point of view 
it is quite natural to focus the interest on the small insta- 
bility branch, as it is precisely anomalous proliferation of 
almost stable segments of trajectories that induces tran- 
sition from exponential to power-law decay of temporal 
correlations. 




FIG. 2: Decay of M^{n) for the system (O with z = 2.0, for 
different cutoffs A, where the red (gray) curve is a regression 
fit of power-law decay with exponent ^ = 1 . 

In fig. ([3]) we show the shape (bimodal) of (normalized) 
finite time distribution P„(A„) for the map as ergod- 
icity predicts, they tend to a Dirac S centered on A in the 
asymptotic limit, but the tails are polynomially "fat" for 
finite times. In the inset (inside fig. ^) we can detect 
an interesting feature: by analyzing the distribution of 
Lyapunov exponents closer to zero in logarithmic scale 
it is possible see how the distribution obeys a power-law 
decay with time. 
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FIG. 3: Distribution of Lyapunov exponents Pn(A„) for the 
map ([6]) with z = 2.0. Inset: a magnification of Pn{\n) is 
shown in logarithmic scale. The distribution of Lyapunov 
exponents closer to zero show a polynomial decay with time. 

As a final remark we point out that, though large devi- 
ation estimates like ([5]) are symmetric with respect to the 



B. Two-dimensional systems: a family of 
area-preserving maps 

The second dynamical system we consider is a family 
of area-preserving maps on the two-torus [— tt, tt)^, de- 
pending upon two parameters s and 7: 



M, 



Vn+i = Vn + f{xn) mod27r, 

Xn+1 = Xn+ Vn+1 mod 2tT , 



where /(a;„) is defined by 



/(x„) = [xn - (1 - e) sin(a;„)]''. 



(8) 



(9) 



A map of this family was introduced in 12311 . and dif- 
ferent features were analyzed in [24,, 2^ l27j : we recall 
a few of the relevant properties. When £ > the map 
is hyperbolic, while for e = the fixed point at (0, 0) 
becomes parabolic: in such a case dynamics close to the 
fixed point depends upon the value of 7, which plays the 
role of an intermittency parameter: correspondingly the 
decay of correlations is exponential in the former case, 
while in the latter a power law is expected, with an ex- 
ponent depending on 7. In p7j the following law was 
proposed 



C(n) 



-3(7+l)/(37-l). 



(10) 



in the case 7 = 1 a rigorous bound is proved [2^ |C(n)| < 
(while (|10p predicts an exponent —3). 
Our numerical experiments encompass both regimes: 
in fig. ([4^) a hyperbolic parameter choice leads to an 
exponential decay rate for A4^(n), while in the intermit- 
tent case a power law is observed, with an exponent in 
agreement with ([TO)) . see fig. (|Dd). In this case the nu- 
merical results were obtained by using 10^ and 10* initial 
conditions uniformly distributed in the phase space, re- 
spectively. 

The case reported in fig. ^p) allows to scrutinize the 
numerical virtues of our approach with alternative meth- 
ods (that were employed in [131): it provides estimates 
as sharp as the analysis of return time statistics (by em- 
ploying however much less initial conditions), while it 
outperforms direct computation of correlations. 



C. Ensemble of modified standard maps 

The last case we present concerns a general issue, 
namely whether universality properties are exhibited for 
Hamiltonian systems with a hierarchical (mixed) phase 
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preserving maps, which we also utilize in the present 
framework. 

The ensemble of modified standard maps on [— 7r,7r)^ 
considered here are given by 

p . f Pn+l = Pn + Ksin{Xn) + K\ 

^'^^ ' \ Xn+1 = Xn+ Pn+1 

where K and K'^ are the nonlinear parameter and mag- 
netic field, respectively. They are uniformly chosen in the 
intervals K G [tt, 1.2 • tt] and K"! € [0.0, 0.4 • tt]. We have 
used 10^ initial conditions taken uniformly in the square 
xq S [—0.1,0.1], po € [—0.1,0.1], well inside the chaotic 
region. The behavior wc show in fig. ([5]) apparently arises 
after a few (~ 10) random realizations of the maps: the 
overall decay rate with an exponent 0.57 ±0.05 is stable, 
as well as the oscillations that may be noticed in fig. ([5]). 
Wether such oscillations are due to finite size effects, cor- 
rections to scaling, or multifractality (ssj , is a feature that 
in our opinion deserves further investigations. Again we 
remark that the present method has an important advan- 
tage that it needs a "small" number of initial conditions 
and realizations of the map to show the asymptotic decay 
of A4^(n) as compared with traditional procedures using 
direct computation of correlation functions or Poincare 
recurrences. 

10° j j 



FIG. 4: Decay Mj^(n) (symbols) together with a regression 
fit (full lines) for the map (|8} with (a) e = 0.5 and 7 = 1 
(exponential rate decay 0.87 ±0.01 ) and (b) e = and 7 = 1 
(polynomial rate decay 3.05 ± 0.05). Both fits were done by 
starting at n = 3000. 

space. This is a much debated issue (we refer to [H, ll^ 
for relevant references over the last twenty years): in par- 
ticular [2^ suggests, on the basis of a Markov tree model 
[30I ] with random scaling factors for transition probabil- 
ities, the existence of an asymptotic algebraic decay for 
the Poincare recurrences. The authors also discuss in 
detail the numerical difficulties in probing such a claim: 
in particular, following ^31i , ,32.] they advocate that an 
average over different Hamiltonian systems significantly 
reduces the extremely long times that are required for 
a single trajectory to sample fine details of the phase 
space structure: this procedure allows them to check for 
an ensemble of area preserving maps their proposal, and 
the polynomial decay exponent of Poincare recurrences 
is estimated as x — 1-57 ± 0.03, which, according to 
[13, [H, [13, [Hi , suggests for correlations the decay law 



The exponent of algebraic decay of Poincare recurrences 
was also numerically investigated in (29j (with an esti- 
mate X — 1-60 ± 0.05), for a different ensemble of area 
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FIG. 5: Distribution of largest Lyapunov exponents for 40 
realizations of the map (|12|l . The black curve is the average 
over all curves and the red (gray) line corresponds to the fit 
of average. The fit decay with an exponent 0.57 ± 0.05. The 
fit was done by starting at n = 500. 



III. CONCLUSIONS 

We have shown how the analysis of large deviations 
for Lyapunov exponents provides a very efficient way 
to investigate correlation and Poincare recurrences de- 
cay for dynamical systems: in particular when dealing 
with weakly chaotic systems, with associated polynomial 
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correlation decay, the power law exponent has been com- 
puted from the distribution function of finite time Lya- 
punov exponents with remarkable accuracy. 

We also remark that in all our tests this procedure 
performs quite well numerically: for instance in our last 
example the estimate has been obtained with a computa- 
tional effort significantly reduced with respect to former 
simulations [28j, i29| . Another interesting point wc plan 
to investigate in the future is to study the decay law of 
correlation function in higher dimensional systems where 
the nonlinearity is not uniformly distributed along differ- 
ent unstable directions fSJl. 
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